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Abstract: We present a theoretical study on the impact of an interface on the FRET 
efficiency of a surface-bound acceptor-donor system. The FRET efficiency can be modified 
by two effects. Firstly, the donor's electromagnetic field at the acceptor's position is changed 
due to the partial reflection of the donor's field. Secondly, both the donor's and the acceptor's 
quantum yield of fluorescence can be changed due to the interface-induced enhancement of 
the radiative emission rate (Purcell effect). Numerical results for a FRET-pair at a glass- water 
interface are given. 
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1. Introduction 

Forster Resonance Energy Transfer or FRET [1] has become one of the most widely used 
fluorescence-spectroscopic tool in biophysical and biomedical research; see the reviews [2-6]. It allows 
for determining inter- and intramolecular distances of a few nanometers with sub-nanometer accuracy, 
and is routinely used for measuring intramolecular conformation and intermolecular interactions. 
In particular with the advent of single-molecule fluorescence spectroscopy, single-molecule FRET 
(smFRET) has become one of the most important toolbox of single-molecule spectroscopists [7,8]. 

Frequently, smFRET experiments are performed with the molecules immobilized on a surface [9,10]. 
This prevents the diffusion of the molecules out of the focal plane and allows for recording long-lasting 
photon bursts from one and the same molecule. Moreover, such measurements can be performed with 
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a Total-Internal Reflection Fluorescence or TIRF microscope for maximizing the signal-to-background 
ratio. However, the proximity of an interface dividing to materials with different refractive indices 
does change the fluorescence emission behavior of fluorescent dyes. In particular, their radiative 
transition rate from the excited to the ground state is changed, which results in a changed fluorescence 
lifetime. Although this effect is small for interfaces between dielectric materials, the present paper 
aims at quantifying the potential impact of such electromagnetic effects on the accuracy of a smFRET 
measurement. To have precise knowledge about such surface-induced effects on the measurable FRET 
efficiency is important in the context of recent efforts to make smFRET measurements as exact and 
quantitative as possible [11-17]. 

In the past, several studies have been published that focus on the strong change of the FRET 
efficiency due to the coupling of the donor and acceptor fluorescence to the local density of states of 
the electromagnetic field as modified by the presence of a cavity [18] or a micro-resonator [19,20]. Here, 
we give a detailed study of this change due to the much weaker influence of an interface dividing two 
materials of different dielectric constants. 



Before considering the problem of the FRET efficiency in the vicinity of an interface, we have to 
first briefly recall the classical theory of FRET, because the details of this derivation will be needed 
for a close understanding of the modifications induced by the presence of an interface that divides to 
half-spaces with different refractive index. In the second subsection, we will then study the change in 
acceptor excitation rate by the presence of an interface. In the third and last subsection, we will consider 
the interface-induced changes in fluorescence quantum yield of both the acceptor and the donor, which 
are important when calculating FRET efficiencies from the ratio of donor and acceptor intensities. 

2.1. Classical Forster Theory 

A light emitting molecule can be fairly well described by the classical picture of an oscillating 
electric dipole. Such a dipole is characterized by the charges q that oscillate, the oscillation amplitude 
and orientation described by the dipole vector a, and the oscillation circular frequency to, which is 
also the frequency of the emitted light. The product of the charge q times the amplitude vector a is 
called the dipole moment (or, more precisely, the amplitude of the dipole moment), i.e., = ga. 
When plugging such an oscillating electric dipole as a source into Maxwell's equations, one obtains 
the electromagnetic field of the well-known Hertzian dipole emitter. In particular, the electric field 
amplitude reads 



2. Theory 
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where r is the distance vector from the dipole to the observation point where the field is calculated, r is 
its modulus, k is the modulus of the wave vector, k = ncu/c, with c being the vacuum speed of light, 
e r is a unit vector pointing from the dipole towards the position r, and n is the refractive index of the 
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medium. Here, and in all that follows, we omit the factor e luJt describing the temporal evolution of the 
electric field. The far-field emission rate of such a dipole is given by 

Sd = \p 2 d kQcn (2) 

where k 0 is the length of the wave vector in vacuum (k 0 = 2n/X), with A being the vacuum wavelength 
of light. When a second molecule (also assumed to be an electric dipole oscillator) is brought very close 
to the emitting molecule (the donor), the field contributions falling off proportional to r _1 and r~ 2 can 
usually be neglected if kr <C 1, and one has to take into account only the near- field component 

E d (r) w [3e r (e r • p d ) - p d ] — — (3) 

n r 

The excitation rate of the second molecule (the acceptor) is given by the product of its absorption cross 
section a a and the number of (virtual) photons per area per second coming from the donor. This photon 
flux is given by the electromagnetic energy flux density associated with the donor's electromagnetic 
field divided by the energy per photon. This energy flux density, averaged over one oscillation period, is 
equal to the modulus of the Poynting vector or (nc/87r) • Ej, and the energy per photon is equal to hc/X, 
where h is Planck's constant. Furthermore, one has to take into account that the acceptor's excitation 
rate is also orientation-dependent and proportional to the scalar product of its dipole vector p a and the 
vector of the exciting electric field E rf . One minor complication comes in by the fact that when one 
measures absorption cross-sections a a on ensembles of molecules, one averages over all possible dipole 
orientations, so that the absorption cross-section of a molecule with its dipole aligned with the electric 
field vector of the exciting field is three times larger than the measured value of a a . Taking all that 
together, one arrives at the following expression for the acceptor's excitation rate k e : 

k e = (3a a )^A| Pa . E(i | 2 (4) 

o7T hC 

where p a is the unit vector along the acceptor's dipole moment. Inserting the explicit expression for 
E d yields 

3a a Xp 2 d 2 3a a Xp d 2 

K = 13 (£V ' Pa) {Gr ' Pd) ~ {Pa ■ Pd)l = 8M^ K (5) 

which also defines the so-called orientation factor k. If both donor and acceptor can freely rotate, one 
can average the orientations p a and p d of donor and acceptor dipole, resulting in k 2 = 2/3. 

The equation for k e still contains the a priori unknown donor dipole amplitude p d . It can be found by 
scrutinizing the total energy emission of the free donor. From a quantum mechanical point of view, the 
totally emitted energy is that of one photon, hc/X. From a classical electromagnetic point of view, it is 
given by the energy emission rate of the donor times the mean radiation-related lifetime of the donor's 
excited state (i.e., its lifetime in the absence of all non-radiative de-excitation channels), which is given 
by the inverse of the radiative transition rate k^ d from the excited to the ground state of the donor. The 
inverse of the sum of the radiative and the non-radiative transition rate, k nr . jd , determines the observable 
excited state lifetime r d , whereas the ratio of kf td and kf jd + k nrjd define the donor's fluorescence 
quantum yield <f) d (i.e., the probability that the donor's energy is given away electromagnetically, and 
not by intermolecular collisions etc.). Thus, kj d is given by the ratio of the observable excited state 
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lifetime r d divided by the fluorescence quantum yield <p d . Setting both the quantum-mechanical and the 
classical-electromagnetic energy emission equal to each other, one obtains 

A3 <p d 

where we have used Equation 2 for the far-field emission rate S d . Thus, one finds 

2 3Hd 



so that the acceptor's excitation rate now reads 



(7) 



_ 9a a <p d 2 

Up to this point, donor emission and acceptor excitation was assumed to happen at exactly one and the 
same wavelength. To obtain the actual acceptor excitation rate, one has to integrate the found expression 
over all wavelengths, taking into account the wavelength-dependence of the donor's emission strength as 
given by its emission spectrum F d (X), and of the acceptor's excitation strength as given by its absorption 
spectrum er a (A). This leads to 

= 9<^ 2 fciAF d (AK(A)A 4 
e 87r(27r) 4 r d n 4 r 6 J d\F d (X) 1 j 

where F d {\)/ J d\F d (\) is the normalized emission spectrum of the donor, and where we have replaced 
k 0 by 2tc/X. The so-called Forster radius is defined by 

90^ 2 f d\F d (\)a a (\)\* 
0 8vr(27r) 4 n 4 J d\F d (X) K ' 



K=-(^) (11) 



so that one finally finds 



Td \ r 

Taking into account the relation between absorption cross section a a (in cm 2 ) and molar extinction 
coefficient e (in cm -1 M _1 ), 

10 3 In 10 

<r a = — ^7 e a (12) 

where Na is the Avogadro-Loschmidt number, one finds the usual textbook expression for R 0 

9OOOlnlO0 d /t 2 /rfAF tj (A)6 Q (A)A 4 
0 128vr 5 iV A n 4 / d\F d {\) U } 

Finally, we have to discuss the two most common methods of experimentally determining FRET 
efficiencies. The first one measures the fluorescence intensity of the donor and the acceptor upon donor 
excitation. For the sake of simplicity we assume that one can neglect direct excitation of the acceptor by 
the excitation light, and any bleed-through of donor emission into the acceptor detection channel. Also, 
the detection efficiencies for both donor and acceptor fluorescence are assumed to be the same. Then the 
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observable fluorescence intensity from the donor in the presence and in the absence of the acceptor will 
be proportional to 

ha 00 k + k k d +fc ° 4) 



and 



Id oc = 0 d , (15) 

respectively, and that of the acceptor is proportional to 

I OC ^f> a ^ 6 ~~ $afce (16) 

^/,a "I" ^nr,a ^/,d "I" ^nr,d "I" k e &/,d "I" ^nr,d "I" k e 

Here, the A;/ and k nr are the radiative and non-radiative transition rates, respectively, and the additional 
subscript d refers to the donor and the subscript a to the acceptor. Accordingly, 0 a is the fluorescence 
quantum yield of the acceptor. Thus, we find the textbook relation for the FRET efficiency E: 

_ Id — ha _ Igj 4>a _ K. _ Rp 

Id I a /4>a + Ida/4>d k f)d + k nr4 + k e .Rfj + r 6 

where we have used Equation 1 1 and the fact that the free donor's de-excitation rate is the inverse of its 
lifetime, k^ d + k nr ^ d = T d x . The second method is based on the measurement of the donor's excited 
state lifetime. In the presence of the acceptor, this lifetime is given by 

Thus, one finds a second relation for the FRET efficiency E as 

£ _ T d — Tda _ ke _ Rp 

Td k ftd + k nr ,d + k e R% + r 6 

2.2. Donor Electric Field above the Surface 

Consider the system depicted in Figure 1: A donor- acceptor pair (smFRET system) are coupled 
together by some rigid ruler keeping them a fixed distance L apart. One of the molecules (let us assume 
the donor) is kept at a distance z d from a surface; the half-space below that surface has the refractive 
index n g , and the half-space above the surface, containing the smFRET system, has refractive index n m . 
As before, both the donor and the acceptor are considered to be ideal electric dipole emitters/absorbers. 
Furthermore, we will assume that both donor and acceptor are free to rapidly rotate around their 
center positions so that one can average over donor and acceptor orientation when calculating FRET 
efficiencies. As will be seen, an extension of the theory beyond this assumption is straightforward but 
would introduce another four degrees of freedom. But we will assume that the connecting line between 
acceptor and donor keeps a fixed angle towards the vertical. For example, in case of a charged surface 
and a charged smFRET system, one of two limiting cases may be expected: either both molecules are 
adsorbed to the interface and 9 = n/2 due to attractive electrostatic interaction, or one finds 9 ~ 0 due 
to electrostatic repulsion between surface and smFRET system. Our results will also be applicable to the 
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more general case of orientational flexibility of the donor-acceptor axis. One has then only to average 
over different orientations of this axis with an appropriate weight function that describes the chance that 
a particular orientation occurs. 

Figure 1. Geometry of the considered problem: A donor molecule with dipole moment 
Pd is kept at a distance z d above a surface. An acceptor molecule with dipole moment p a 
is situated a fixed distance L away from the donor, the angle between the vertical and the 
connecting line is 9. 




x 



As we have seen in the preceding subsection, the energy transfer rate k e is proportional to the 
donor-generated local electric field strength at the position of the acceptor. Due to partial reflection at 
the interface, this field strength will be different from the situation without interface. For calculating this 
field strength, we start by expanding the electric field of the free donor into a Weyl representation [21] 
of plane waves (see Figures 1 and 2), which is the standard field representation when considering 
planar problems, similar to using Bessel function expansions for problems with cylindrical symmetry, or 
spherical harmonics expansions for problems with spherical symmetry: 

MP, z) = ^ jj ^ [ef (ej ■ p d ) + e s (e s ■ p d )] e ^+H*- rf | (20 ) 

where q is the part of the wave vector parallel to the interface, and w = ^k 2 — q 2 its part perpendicular 
to the interface. The unit vectors and e s are pointing along the electric field polarization of p- and 
s-waves, respectively, as shown in Figure 2, where the plus sign applies for z > Zd, and the minus sign 
for z < z d . The integration over q extends over the whole two-dimensional q-plane. The electric field 
as given in Equation 20 is the field of an oscillating electric dipole in free space with refractive index n m 
and yields, after explicit integration, the result of Equation 1. 
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In the presence of an interface, each plane wave component in the Weyl representation is partially 
reflected, adding to the total field strength in the upper half space. Using the standard, g-dependent 
Fresnel coefficients r P;S for the reflection of plane p- and s-waves at a planar interface, the reflected field 
can be written as 

~2 r r A~ r n 

e iw{z+z d )+i<\-p ^21) 



E, 



P, z ) = / / — [ r p(<l) e t ( e P ■ Pd) + r s(q) e s (e, • p d ) 



so that the total electric field strength in the upper half space is given by E' rf = + E r . Whereas 
the integrals in Equation 20 can be calculated analytically, the integrals in Equation 21 have to be 
calculated numerically. The excitation efficiency of the acceptor by the field of the emitting donor is 
then proportional to the square of the scalar product between electric field amplitude and acceptor dipole 
moment, i.e., proportional to \E' d (p a , z a ) ■ p a | 2 , where p a and z a are the lateral and vertical position of 
the acceptor, respectively, and p a its excitation dipole moment. What we are interested in is a comparison 
of the FRET efficiency of our acceptor-donor system in free space with the same system attached to the 
surface. For the sake of simplicity, we will here compare both situations in the limit of rapidly rotating 
acceptor and donor orientation. As already mentioned, we could equally well consider particular fixed 
orientations of donor and acceptor, which would, however, add another four degrees of freedom without 
making the general effect of an interface on the FRET efficiency more clear. As will be seen from 
the calculations of the emission rate enhancement below, all the considered electrodynamic effects are 
strongly orientation-dependent, so that a fixed donor and/or acceptor orientation will lead, in general, to 
an even stronger modification of the FRET efficiency than in case of rapid rotational diffusion of donor 
and acceptor. Thus, we define the ratio 

Q(P«,* a )= { ™ Pa,Z i' P i\ n (22) 

(\Ed[P a ,Za) -Pa\ )fl 

which tells us how much the acceptor's excitation rate of the surface-bound system is enhanced (Q > 1) 
or reduced (Q < 1) when compared with a measurement in free solution. Here, the angular brackets 
with subscript indicate averaging over all orientations of donor and acceptor. Due to symmetry, the 
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ratio Q depends only on the modulus p a = \p a \ of the lateral position vector, but not its two-dimensional 
orientation. Moreover, the problem is reciprocal: the ratio Q in Equation 22 remains the same when 
interchanging the positions of donor and acceptor. This can be directly checked by repeating the above 
calculations with interchanged positions. 

2.3. Emission Rate above the Surface 

Besides a change in the excitation intensity felt by the acceptor due to the back-reflection of the 
donor's electromagnetic field, there is a second surface-induced aspect affecting the FRET efficiency. As 
we have seen in Section 2.1, the acceptor excitation rate k e is also proportional to the donor's radiative 
transition rate kf d = 4> d / r d- During the emission of the donor, the back-reflected electromagnetic field 
acts on the emitting donor itself, enhancing or inhibiting its radiation and thus changing its radiative 
transition rate. This is also called the Purcell effect after Edward Mills Purcell who first pointed it out 
in 1946 [22]. When the reflected field E r acts on the oscillating dipole, the dipole performs the average 
work per time of (see [23]) 

W = ^ Pd -Q(E r ) (23) 

which adds an additional term to the energy loss per time exhibited by the dipole. Thus, the emission 
rate is changed to 

S' d = S d + W (24) 
so that the relative change 7 in emission rate alias radiative transition rate is given by 

^trf^ i+ ^<^ (E ^ (25) 

where we have again averaged over all possible orientations of the emission dipole, p^. The factor 7^ is 
the interface-induced change of the radiative transition rate of the donor. An identical relation holds also 
for the emission rate of the acceptor, where one has only to exchange all subscripts d by a. If one can 
neglect the optical dispersion of the materials in the upper and lower half spaces, then the resulting factor 
7 is identical for both donor and acceptor if they are at the same distance from the surface. Knowing 
both the factor Q and 7, the interface-modified acceptor excitation rate is given by 

K = ~?Qk e . (26) 

The modified radiative emission rate changes also the quantum yield of both the donor and the acceptor, 
according to the formula 

_ k 'f ikf 70 



k' f + k nr "/kf + k nr 70 + 1-0 
where we have omitted the subscript d or a. 



(27) 
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2.4. FRET Efficiency above the Surface 

When calculating the FRET efficiency E by the intensity ratio according to Equation 17, two scenarios 
are possible. If one uses the interface-modified values of the fluorescence quantum yields of both donor 
and acceptor, then the interface-modified FRET efficiency is given by 

E , = Wa = 1_QE 

la/Pa + IM 7Q£+(70d + l-0d)(l-£) 



However, if one uses the quantum yield values for the dyes in free solution, the calculated FRET 
efficiency is 

, _ LI (Pa _ iQE 

E = hi + Ida/fa = 7Q£+(70a + l-0a)(l-£) ^ 



which becomes identical to Equation 28 only if <p a = <p d . Similarly, when using the lifetime-based 
Equation 19 for determining the FRET efficiency, one can either use the lifetime value r<2 of the donor 
in free solution, or the interface-modified value r' d = r&j {^4>d + 1 — <j>d)- When using r' d , one recovers 
Equation 28, and when using r d , one finds 

E , = iQE + ^-l^-E) 

7 QE+(j<p d +l-<j> d )(l-E) K ) 



3. Results and Discussion 

We first performed calculations of 7, Q, and E' for a glass-water interface, where the refractive index 
of glass was set equal to 1.52 and that of water equal to 1.33. All calculations where done for an assumed 
emission wavelength of 550 nm. Figure 3 shows the distance dependence of the radiation enhancement, 
7. As can be seen, close to the interface (a few dozen nanometers), 7 does not change significantly, 
reaching a maximum value of 1.18 on the surface (averaged over all possible dipole orientations). For 
the sake of completeness, the figure shows also the extreme values of 7 for horizontal (||) and vertical 
(_!_) emission dipole orientation. An identical figure is obtained for the radiation enhancement of an 
acceptor, where the x-axis has only to be rescaled by the ratio of the acceptor's to the donor's emission 
wavelength. Taking into account the dependence of the fluorescence quantum yield on the enhancement 
factor 7, Equation 27, we find for the relative change in quantum yield A0/0 = (1 — 0)A7/7, so that a 
18% change in the radiative rate 7 for a fluorophore with 50% fluorescence quantum yield translates into 
an approximately 9% change in quantum yield. Of course, the effect becomes smaller with increasing 
value of quantum yield. 
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Next, Figure 4 presents the dependence of both the factor Q(p a , z a ) and the FRET efficiency E' as a 
function of lateral (p a ) and vertical (z a ) distance between donor and acceptor, respectively, where it is 
assumed that the donor sits directly on the surface at pd = 0. As was already mentioned, an identical 
figure is obtained when exchanging the positions of donor and acceptor. 

The Q-factor varies between 0.87 and 1.06, depending on the orientation of the donor-to-acceptor 
axis. For example, for a donor- acceptor distance of 10 nm, if the donor-acceptor connecting line is close 
to the surface, Q is close to its minimum of 0.87, whereas if this axis is perpendicular to the surface, it 
reaches its maximum value of 1.06. When using Equation 28 or 29 for calculating the FRET efficiency, 
then the relative change in E is given by 



^ = (1 _ B) ^ + (1 _ B)(1 _^ (3!) 

where <\> = <\> d for Equation 28, and <\> = <\> a for Equation 29. Thus, at midpoint FRET efficiency of 0.5, 
the relative change AE / E is half that of AQ (Q = 1) for dyes with quantum yield of 1.0, and the FRET 
efficiency is only changed between —7% and +3%. For quantum yield values of 0.5, the relative change 
in FRET efficiency lies between —2.5% and 4.5%. Thus, although the energy transfer rate k e can vary 
by up to 13%, and the radiative emission rate by up to 18%, the resulting variation in FRET efficiency is 
much smaller due to the intricate dependence of E on both these factors. 
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Figure 4. Enhancement Q of the energy transfer rate (left) and the FRET efficiency (right) 
as a function of the relative position of the acceptor (donor) if the donor (acceptor) is 
positioned directly on a glass surface (n = 1.52). For the FRET efficiency calculations 
(right panel), it was assumed that its value in free solution is 0.5 and that the quantum yield 
of the fluorophores is 1.0. 
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In recent years, objectives with N.A.-values of 1.65-1.7 have become available and are used for 
single-molecule detection and spectroscopy experiments, because they provide exceptional high light 
collection efficiency and enable efficient TIR illumination with very small penetration lengths of the 
evanescent field. When using these objectives, one needs special sapphire cover slides and high-index 
immersion oils with a refractive index of 1.78. We repeated all the above calculations for checking how 
such a high refractive-index material will influence the smFRET measurement. For a fluorophore close 
to the sapphire/water interface, the radiative emission rate is now enhanced by 45% (averaged over all 
dipole orientations) in contrast to the 18% for the glass/water interface. Figure 5 shows the results for 
both the transfer rate Q and the FRET efficiency E' for the same geometry and distances as used in the 
calculations for the glass/water interface. The enhancement factor Q for the transition rate varies now 
between 0.71 and 1.2, and the resulting FRET efficiency varies between —17% and +9% for dyes with 
quantum yield 1.0, and between —9% and +17% for dyes with quantum yield 0.5. Thus, although the 
impact of a glass/water interface on the FRET efficiency is rather moderate, it is no longer negligible for 
a sapphire/water interface. 

There are two possible and straightforward ways to experimentally test the predictions made in this 
paper. The simplest way would be to compare FRET measurements of the same system in free solution, 
and when it is bound to a surface. However, a problem might be that other factors like the different 
nano-environment close to a surface or the binding chemistry used for surface attachment may also 
influence the measured FRET efficiency. A more elegant way would be to use a short and rigid DNA 
double strand with donor and acceptor at both ends as the smFRET system, and to bind it with one end to 
the surface. If the surface is prepared in such a way that it carries a negative net charge, one can assume 
that the negatively charged DNA double strand will adopt an orientation orthogonal and away from the 
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surface, whereas for a positively charged surface, the double strand will adsorb to the surface and thus 
lie flat on it. As can be seen from Figures 4 and 5, one expects an opposite effect on the FRET efficiency 
for both cases: For an orientation perpendicular to the surface, the FRET efficiency is enhanced, whereas 
for an orientation parallel to the surface, it is decreased. At least for a sapphire substrate, the difference 
between both situations (±10%) should be easily observable. 

Figure 5. Enhancement Q of the energy transfer rate (left) and the FRET efficiency (right) as 
a function of the relative position of the acceptor (donor) if the donor (acceptor) is positioned 
directly on the sapphire surface (n = 1.78). For the FRET efficiency calculations (right 
panel), it was assumed that its value in free solution is 0.5 and that the quantum yield of the 
fluorophores is 1 .0. 
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4. Conclusions 

In the present paper, we have given a complete theoretical analysis of the impact of interface-induced 
changes in the donor's electric field as well as the radiative transition rates of both donor and acceptor 
on the FRET efficiency. Although the donor's electric field as well as the radiative transition rates can 
change by more than ten percent, the resulting impact on FRET efficiency for a glass/water interface 
is rather moderate. This is, however, no longer true when measurements are performed on high 
refractive-index substrates such as sapphire, which is used in some single-molecule experiments. Finally, 
we have given some indications of how to experimentally measure the predicted effects. 
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